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Abs t rac t 

Droplet size and velocity and gas velocity 
were measured in a water spray using a two- 
component Phase/Doppler Particle Analyzer. A com- 
plete set of measurements was obtained at axial 
locations from 5 to 50 cm downstream of the noz- 
zle. The nozzle used was a simple axisymmetric 
air-assist nozzle. The sprays produced, using the 
atomizer, were extremely fine. Sauter mean diame- 
ters were less than 20 pm at all locations. Meas- 
urements were obtained for droplets ranging from 1 
to 50 pm. The gas phase was seeded with micron 
sized droplets, and droplets having diameters of 
1 4 pm and less were used to represent gas-phase 
properties. Measurements were compared with predic- 
tions from a multi-phase computer model. Initial 
conditions for the model were taken from measure- 
ments at 5 cm downstream. Predictions for both the 
gas phase and the droplets showed relatively good 
agreement with the measurements. 

Nomenclature 

a acceleration of gravity 

C[ constants in turbulence model 

Cd drag coefficient 

d nozzle diameter 

dp droplet diameter 

k turbulence kinetic energy 

L e dissipation length scale of eddy 

nip droplet mass 

a number of droplet groups 

ni number of droplets per unit time in group i 

Re Reynolds number 

r radial distance 
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S source term 

Sp source term due to droplets 

t time 

t e eddy lifetime 

u axial velocity 

Vj volume of computational cell j 

v radial velocity 

w angular velocity 

x axial distance 

Xp droplet position vector 

e rate of dissipation of turbulence kinetic 

energy 

p viscosity 

Pt turbulent viscosity 

p density 

i turbulent Prandt 1 /Schmidt number 

4 j generic property 

Vxp relative path length of droplets in an eddy 
7tp droplet residence time in an eddy 
Subsc ripts : 

c centerline quantity 

o initial axial location 

p droplet property 

® ambient condition 

Superscripts : 

( ) fluctuating quantity 

(~) time averaged value 

(+) vector quantity 

(~) Favre averaged value 
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Introduction 


Spray nozzles have many important applications 
and have been the subject of considerable research. 
The objective of the present investigation was to 
obtain measurements of droplet size, velocity, and 
volume flux and gas velocity in the flowfield of a 
simple atomizer and compare their values with pre- 
dictions from a multiphase computer model. 

A similar investigation was performed by Dodge 
and Schwalb 1 using a pressure swirl atomizer. Most 
of the initial conditions for their model were esti- 
mated, however, predictions were shown not to be 
very sensitive to initial conditions. For that 
study, air velocities were very low in comparison 
to droplet velocities. A single-component Phase/ 
Doppler instrument was used for the measurements. 

The liquid utilized was an aircraft fuel system cal- 
ibration fluid (MIL-C-7024 type II). The computer 
model used for comparison with the data was the 
FLUENT code. In an earlier study, reported by 
Solomon, et a 1 . , 2 an air-atomizing nozzle operating 
with nonevaporating vacuum pump oil was investi- 
gated. Gas-phase velocities were measured with a 
standard LDA system. Droplet size and velocity 
were measured using a double-flash shadow photo- 
graphic method. Droplet fluxes were measured using 
an isokinetic sampling probe. Experimental measure- 
ments were compared with predictions from three two- 
phase flow models . 


Other experimental spray studies using a 
Phase/Doppler instrument are reported in Refs. 3 
and 4. Measurements are presented for both droplet 
and gas-phase velocities. Presser and Semerjian 5 
also report a spray study using an ensemble light 
scattering technique. No gas-phase measurements 
are presented. A recent study is reported in Ref. 6 
describing Phase/Doppler measurements of a water 
spray in a swirling flowfield. Both gas phase and 
droplet measurements are reported. 


Experimental Methods 

The sprays were directed vertically downward 
within a large enclosure (1.8 by 1.8 by 2.4 m high). 
The air-assist nozzle could be traversed vertically 
within the enclosure while the enclosure and nozzle 
could be traversed together in two horizontal direc- 
tions. This arrangement allowed rigid mounting 
of all optical instrumentation used during the 
study. The flow downstream impinged upon a screen 
located 1 m below the measurement plane to remove 
larger droplets and allow smaller droplets to be 
used as seeding for the ambient air. 

The air-assist nozzle is shown in Fig. 1. The 
inside diameter of the liquid tube was 0.39 mm and 
the orifice for the nozzle was 3.18 mm. This noz- 
zle was selected because it is relatively simple 
and does not utilize a swirler for atomization. 

Water was used for the liquid spray and was supplied 
using a tank pressurized with nitrogen to provide a 
steady flowrate without pressure oscillations. Liq- 
uid flowrates were measured using a mass flowmeter. 
Air flowrates were measured using a 1.27 mm diame- 
ter cri t ical-f low orifice. The air flowrate was 
controlled by regulating upstream pressure of the 
air. The atomizing air was seeded with micron 
sized droplets supplied from a Norgren lubricator. 
The fluid used for the seeding was Rosco fog fluid, 
typically used for smoke generators. The average 


measured drop size in the jet using only the seeder 
was 1.5 pm. For the results reported in the present 
paper, the air and water flowrates were 2.25 and 
0.49 g/s, respectively. 

All measurements were performed using a two- 
component Phase/Doppler Particle Analyzer. Previous 
studies describing this instrument have been 
reported in Refs. 7 and 8. The configuration cho- 
sen for the present experiment utilized a 2-W argon- 
ion laser, frequency shifting for both channels 
using rotating diffraction gratings, a 500 mm focal 
length focusing lens, and a 450 mra focal length 
receiving lens. Light was collected 30° off-axis 
in the forward direction and focused onto a slit. 

For measurements at 5 and 10 cm downstream of the 
nozzle, a 50 pm wide slit was used. The remainder 
of the measurements used a 100 pm wide slit. Gener- 
ally, 64 000 total samples were taken at each meas- 
urement location. The number of validated samples 
varied from about 90 to 20 percent of the total, 
depending on the region of the flowfield where meas- 
urements were being obtained. Measurements in loca- 
tions near the center of the spray close to the 
injector exit where number densities were extremely 
high and velocities were also high showed the larg- 
est number of signal rejections. Regions downstream 
of the nozzle where the spray was not as dense and 
velocities lower had much higher validation rates. 
Droplet diameters ranging from 1 to 50 pm were meas- 
ured. Velocity measurements are reported for drop- 
let diameters of 4, 11, 18, and 25 pm. Gas-phase 
velocities are obtained from measurements of drop- 
lets with diameters of 1.4 pm or less, A traverse 
across the entire spray was performed at each axial 
measurement station. 


Theoretical Methods 

The analysis is limited to steady, ax i symmet- 
ric, dilute, nonevaporating, droplet- laden , non- 
swirling, turbulent jets in an infinite stagnant 
medium. The boundary- layer approximations are 
adopted in conjunction with a k-c turbulence 
model for closure. The assumptions for the continu- 
ous phase are identical to those reported in Ref. 9. 
Favre-averaged equations (equivalent to time- 
averaged in this case) govern the variation of mean 
axial and radial velocities, turbulent kinetic 
energy, and rate of dissipation of turbulent kinetic 
energy. With these assumptions, the governing equa- 
tions for the gas phase can be expressed in the fol- 
lowing general form: 


f^(pu<t>) + p §prpv<t>) 
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The parameters and S* appearing in Eq. (1), 
as well as empirical constants used, are listed in 
Table I. 

The turbulent viscosity was calculated as: 
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Boundary conditions for Eq. (1) are: 
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Initial conditions were obtained from the meas- 
urements obtained at 5 cm downstream of the nozzle. 
c 0 was calculated From the definition of a turbu- 
lent length scale and is given as follows: 
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where, nj is the number of particles per unit time 
in each group. The calculations for the continuous 
phase were performed using a modified version of 
Genmix. 11 The computational grid utilized has 
33 cross-stream grid nodes and streamwise step size 
was limited to 2 percent of the current flow width 
or an entrainment increase of 2 percent - whichever 
was smaller. For the droplet calculations, 19 200 
droplet groups were tracked through the flowfield. 

Results and Discussion 


where L was taken as 8 percent of the flow width 
at that axial location. Initial values of w’ 
were assumed to be equal to measured values of v' . 


The computation of the liquid phase is based 
on the solution of Lagrangian equations of motion 
for a statistically significant number of repre- 
sentative droplet groups. Major assumptions for 
the droplet trajectory calculations are: quasi- 

steady gas phase, no droplet-droplet interaction, 
and use of a solid-sphere drag law. With these 
assumptions, the position and velocity of each drop- 
let group can be found by integrating: 
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where 1*1.2, 3 for each one of the Cartesian 
coordinates and the velocities shown in these equa- 
tions are instantaneous velocities. The drag coef- 
ficient was calculated as follows: 
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c d - sK l * ! V). Re < 1000 ' 

C D - 0.44, Re > 1000 (7) 

The effects of turbulent gas phase fluctuations 
on droplet trajectory are incorporated using a sto- 
chastic separated-f low (SSF) model. The interaction 
between turbulence and droplets is considered by 
using a random-walk technique. Each group of drop- 
lets interacts with a succession of turbulent 
eddies. The interaction time is taken to be the 
minimum of either the eddy lifetime or the transit 
time for the droplet to cross the eddy. Properties 
of each eddy are found by making a random selection 
from the pdf of velocity. Instantaneous eddy veloc- 
ity is calculated assuming velocity fluctuations 
are isotropic with a Gaussian pdf having a most 
probable value found from the local mean velocity 
and a variance of 2k/3 . Details of the procedure 
can be found in Refs. 2 and 10. The interaction 
between the droplets and the gas phase yields 
source terras in the governing equations for conser- 
vation of momentum. They are found by computing the 
net change in momentum as each droplet group i 
passes through a computational cell j and is given 
as follows: 
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Axial Profiles 


The variation of centerline values with axial 
distance are illustrated in Fig. 2(a) to (e). Both 
measurements and predictions are presented for axi- 
al velocities of both gas and liquid phases. As 
illustrated in Fig. 2(a) to (d) , axial velocities 
at 5 cm downstream of the nozzle are fairly high, 
ranging from 80 m/s for the gas phase to 93 m/s for 
the 25 pm drops. As measured, the droplets had 
higher velocities than the gas at this axial loca- 
tion. Axial velocities decay very quickly from 5 
to 20 cm downstream of the nozzle for both the gas 
phase and the three droplet sizes shown. Results 
for the 4 pm droplets were almost identical to the 
gas phase and are not shown. Even though the drop- 
let size is relatively small, larger droplets pos- 
sess higher axial velocities at distances less than 
20 cm from the nozzle. At distances greater than 
about 20 cm, velocity differences between droplet 
sizes are considerably reduced. As illustrated in 
Fig. 2, predictions from the model for axial center- 
line velocities are in reasonably good agreement 
with the measurements. At downstream distances of 
5 and 10 cm, the centerline measurements presented 
for the 25 pm droplets are based on less than 50 
samples, consequently they should be used with cau- 
tion. Sample sizes for the 18 pm and smaller drops 
were adequate at all downstream locations. At lar- 
ger downstream distances, centerline measurements 
for the 25 pm droplets are based on at least 2000 
samples. Figure 2(e) presents the axial centerline 
profile of turbulence kinetic energy for the gas 
phase, normalized by the local centerline axial 
velocity. For the experimental values, since only 
u’ and v f were measured, it was assumed that w’ 
was equal to v f . The increase of normalized turbu- 
lent kinetic energy with axial distance is initially 
overpredicted, however, predictions show better 
agreement with the measurements at increasing axial 
distances from the injector. Mo attempt was made 
to account for the modification of the gas-phase 
turbulence by the droplets. 

Radial Profiles 


Radial profiles of spray properties at 10 cm 
downstream are illustrated in Fig. 3(a) to (f). 
Results are presented for the 18 pm diameter drop- 
lets and the continuous phase. Radial profiles for 
the other droplet sizes were similar except for the 
fluctuating quantities, which will be discussed 
later. Gas phase r adia l profiles of mean axial 
velocity, k, and u'v' are shown in Fig. 3(a) to 
(c). Centerline axial velocity has decayed to 
about 43 m/s and is underpredicted by abo ut 2 0 per- 
cent. Centerline k and peak values of u’v are 


3 


somewhat overpredicted at this axial location, how- 
ever the jet width is fairly well predicted. Pre- 
dicted and measured radial profiles of mean axial 
velocity, and fluctuating axial and radial veloci- 
ties of 18 pm diameter droplets are illustrated in 
Fig. 3(d) to (f). Again, the peak axial velocity 
is underpredicted, but only by about 10 percent. 
Fluctuating axial velocities are substantially 
greater than fluctuating radial velocities at this 
axial location. Predicted values of fluctuating 
axial velocity show lower values than the measure- 
ments, while fluctuating radial values agree reason- 
ably well with the measurements. This is due to 
the assumption of isotropic fluctuating velocities 
in the model, while measured axial fluctuating 
velocities were larger than radial fluctuating 
velocities. The predicted values show some scatter 
in the profiles. Predictions for the 18 pm droplets 
were obtained using 4800 droplet groups. Increasing 
the number of droplet groups would result in 
smoother radial profiles. Larger droplets showed 
larger fluctuating axial velocities than smaller 
droplets near the centerline at this axial location, 
however, this may be due to the fact that fewer 
larger droplets are found at the center of the spray 
and the sample size may not be large enough. An 
examination of the complete data set showed that 
larger droplets consistently showed lower fluctuat- 
ing radial velocities across the entire width of 
the spray than smaller droplets. This was also 
shown by the predictions. Fluctuating radial veloc- 
ities are important for turbulent dispersion of the 
droplets since they are larger than the mean radial 
velocities in the present spray. 

Radial profiles at 20 cm downstream of the noz- 
zle are illustrated in Fig. 4(a) to (f). Results 
are again presented for the gas phase and 18 pm 
diameter droplets. Figure 4(a) to (c) presents gas 
phase measurements and predictions of mean axial 
velocity, k, and u’v'. The centerline gas veloci- 
ty has decayed to about 22 m/s from the value of 
80 m/s measured at 5 cm downstream of the nozzle. 

The predictions agree reasonably well with the meas- 
urements; the centerline is underpredicted by abo ut 
6 percent. Radial profiles of k and also u’v’ 
show fairly good agreement with the measurements at 
this axial location, although k is underpredicted. 
Radial profiles of the 18 pm droplets are presented 
in Fig. 4(d) to (f). Even at 20 cm downstream of 
the nozzle, mean axial velocities of the droplets 
are slightly, but consistently, larger than the gas 
phase. Peak centerline velocity is underpredicted 
by about 4 percent, however values elsewhere show 
very good agreement with the measurements. Fluctu- 
ating axial and radial velocities of the 18 pm drop- 
lets again show effects of anisotropic velocity 
fluctuations, but are reduced here compared to val- 
ues observed at 10 cm downstream of the nozzle. 
Predicted radial fluctuating velocities agree with 
the measurements, while axial fluctuating veloci- 
ties are still lower than the measurements. 

Profiles at 30 cm downstream are shown in 
Fig. 5(a) to (f) for both the gas phase and the 18 
pm droplets. Cas phase axial centerline velocity 
has decayed to about 14 m/s. Predictions of mean 
axial velocity show excellent agreement with the 
meas urements at this location. Gas phase k and 
u’v’ are als o shown, where some scatter is evident 
in the u’v' data. The reason for the scatter is 
probably due to a measurement inaccuracy at this 
location. Droplet profiles at this axial location 
are also presented for the 18 pm droplets. Predic- 


tions for mean axial velocity and fluctuating 
radial velocity show quite good agreement with the 
measurements. Predictions of fluctuating axial 
velocity again are lower than the measurements, how- 
ever, the differences are decreasing as axial dis- 
tance increases and axial velocity decreases. 

Radial profiles at the largest axial downstream 
measurement location of 50 cm are presented in 
Fig. 6(a) to (f). As illustrated in Fig. 6(a) 
to (c), predictions of gas phase radial properties 
show fairly good agreement with the measurements at 
this location, although k is still underpredicted, 
especially near the centerline. The data show a 
slight shift to the positive r direction of about 
2 mm at this axial location. This is probably due 
to a slight misal 1 ignment between the nozzle axis 
and the gravity vector which is more evident at 
larger axial distances. Radial profiles for the 18 
pm diameter drops are presented in Fig. 6(d) to 
(f). Centerline mean axial droplet velocity has 
decayed to slightly less than 10 m/s at this axial 
location. Fluctuating radial velocity is still 
smaller than fluctuating axial velocity, however, 
the flow becomes more isotropic. Predicted values 
show the same distribution trends as the measure- 
ments, but values are on the order of 20 percent 
lower than the measurements. 

Radial profiles for the 4, 11, and 25 pm 
droplets are not presented. Results for the 4 pm 
droplets are almost identical to the gas phase. 

The results for the 11 pm droplets are somewhere 
between the gas phase and the 18 pm droplets. The 
25 pm droplets show larger differences between the 
drops and the gas phase, but there are relatively 
few of these in the spray, especially near the cen- 
ter of the jet near the nozzle. 

Conclusions 

A complete set of measurements and predictions 
of both the gas phase and droplets have been pre- 
sented for a water spray spraying downward into 
stagnanant surroundings. The flow geometry was 
kept relatively simple so that an axisyinmetric , par- 
abolic, computer model could be used. Major conclu- 
sions from the study are listed as follows: 

1. Axial centerline velocity decay with dis- 
tance downstream of the nozzle was reasonably well 
predicted for both the gas phase and the droplets. 
For the 18 pm diameter droplets, droplet mean axial 
velocities decayed from 86 m/s at 5 cm downstream 
of the nozzle to 9 m/s at 50 cm downstream. The 
gas phase decayed from 80 m/s at 5 cm to 9 m/s at 
50 cm downstream of the nozzle. 

2. Predicted radial profiles of mean properties 
for both the gas phase and the droplets showed rela- 
tively good agreement with the measurements at all 
axial stations. Better agreement was obtained as 
axial distance increased. 

3. For the 18 pm droplets, fluctuating axial 
velocities were underpredicted while fluctuating 
radial velocities generally showed good agreement 
with the measurements. This was probably caused by 
the isotropic assumption of the SSF model. 

4. The spray produced by the nozzle was very 
axi symmetric making it useful for comparison with 
model predictions. 
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Figure 1 . — Diagram of air-assist nozzle. 
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(a) Gas phase mean axial velocity. 


(b) Gas phase turbulence kinetic energy. 


(c) Gas phase shear stress. 


(d) Mean axial velocity of 18 micron drops. 
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Figure 6 —Radial profiles at 50 cm downstream of the nozzle. 
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